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Multistatic Cloud Radar Systems: Joint Sensing and 

Communication Design 

Seongah Jeong, Shahrouz Khalili, Osvaldo Simeone, Alexander Haimovich and Joonhyuk Kang 


Abstract —In a multistatic cloud radar system, receive sensors 
measure signals sent by a transmit element and reflected from 
a target and possibly clutter, in tbe presence of Interference and 
noise. Tbe receive sensors communicate over non-ideal backbaul 
links with a fusion center, or cloud processor, where the presence 
or absence of the target is determined. The backhaul architecture 
can be characterized either by an orthogonal-access channel or by 
a non-orthogonal multiple-access channel. Two backhaul trans¬ 
mission strategies are considered, namely compress-and-forward 
(CF), which is well suited for the orthogonal-access backhaul, 
and amplify-and-forward (AF), which leverages the superposi¬ 
tion property of the non-orthogonal multiple-access channel. In 
this paper, the joint optimization of the sensing and backhaul 
communication functions of the cloud radar system is studied. 
Specifically, the transmitted waveform is jointly optimized with 
backhaul quantization in the case of CF backhaul transmission 
and with the amplifying gains of the sensors for the AF backhaul 
strategy. In both cases, the information-theoretic criterion of the 
Bhattacharyya distance is adopted as a metric for the detection 
performance. Algorithmic solutions based on successive convex 
approximation are developed under different assumptions on 
the available channel state information (CSI). Numerical results 
demonstrate that the proposed schemes outperform conventional 
solutions that perform separate optimizations of the waveform 
and backhaul operation, as well as the standard distributed 
detection approach. 

Index Terms —Cloud Radio Access Networks (C-RANs), quan¬ 
tization, localization, Cramer-Rao bound (CRB). 

I. Introduction 

This papeiQ addresses a distributed radar system that in¬ 
volves sensing and communication: a transmit element illumi¬ 
nates an area of interest, in which a target may be present, and 
the signals returned from the target are observed by sensors. 
The sensors have a minimal processing capabilities, but com¬ 
municate over a backhaul network with a processing center, 
referred to henceforth as a fusion center, where target detection 
takes place (see Fig. [T]). Such architecture is different from a 
classical multistatic radar system in which each constituent 
radar performs the full array of radar functions, including 
target detection and tracking. The considered architecture is 
motivated by the proliferation of low-cost, mobile or fixed 
sensors in the “Internet of Things,” which are supported by 
global synchronization services such as the global positioning 
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system (GPS), and are capable of communicating with a fusion 
center in the “cloud” through a backhaul wireless or wired 
network. For example, the receive sensors could be mounted 
on light poles, trucks or unmanned aerial vehicles (UAV’s) 
and could be connected to a wireless access point via Wi¬ 
Fi or dedicated mmWave links. As this architecture can be 
implemented by means of cloud computing technology, we 
refer to it as “cloud radar.” 

The main purpose of this work is to study the interaction 
between the sensing and backhaul communication functions 
in a cloud radar architecture, and to develop an understanding 
of the performance gain to be expected by means of a joint 
optimization of these two functions, namely of waveform 
design for sensing and of backhaul transmission. 

A. Background 

The separate design of radar waveforms, under the as¬ 
sumption of an ideal backhaul, has long been a problem 
of great interest 0, 0. For monostatic radar systems, i.e., 
radars with single transmit and receive elements, optimal 
waveforms for detection in the Neyman-Pearson sense were 
studied in HI. In a multistatic radar system, where the signals 
received by a set of distributed sensors are processed jointly, 
the performance of the Neyman-Pearson optimal detector is 
in general too complex to be suitable as a design metric. 
As a result, various information-theoretic criteria such as the 
Bhattacharyya distance, the Kullback Leibler divergence, the 
J-divergence and the mutual information, which can be shown 
to provide various bounds to the probability of error (missed 
detection, false alarm and Bayesian risk), have been considered 
as alternative design metrics S-Q. 

Instead, the separate design of backhaul communication 
functions, for fixed radar waveforms was studied in 0- 
ca under a compress-and-forward (CF) strategy, for which 
backhaul quantization was optimized, and in ifT^ . ifTSll under 
an amplify-and-forward (AF) scheme, for which the power 
allocation at the sensors was investigated using the minimum 
mean square error (MMSE) as the performance criterion. 

B. Main Contributions 

Unlike prior work, in this paper, we tackle the problem 
of jointly designing the waveform, or code vector, and the 
transmission of the receive sensors over the backhaul. This ap¬ 
proach is motivated by the strong interplay between waveform 
and backhaul transmission designs. For instance, waveform 
design may allocate more power at frequencies that are less 
affected on average by clutter and interference, while the 
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Fig. 1. Illustration of a multistatic cloud radar system, which consists of 
a transmit element, N receive sensors, and a fusion center. All the nodes 
are configured with a single antenna. The receive sensors are connected to 
the fusion center via orthogonal-access or non-orthogonal multiple-access 
backhaul links. 


backhaul transmission strategies are adapted accordingly to 
devote most backhaul resources, namely capacity or power, to 
the transmission of such frequencies to the fusion center. 

Two basic types of backhaul links between the radar receive 
sensors and the fusion center are considered, namely orthog¬ 
onal and non-orthogonal access backhaul. In the former, no 
interference exists between the sensors, as in a wired backhaul, 
while in the latter, the backhaul forms a multiple-access 
channel, where channels are subject to mutual interference, as 
in a wireless backhaul. Furthermore, two standard backhaul 
transmission schemes are investigated, namely CF and AF. 
As in the Cloud Radio Access Network (C-RAN) architecture 
in communication m, CF is particularly well suited to an 
orthogonal backhaul architecture: each sensor satisfies the 
backhaul capacity constraint quantizing the received baseband 
signals prior to transmission to the fusion center. AF, instead, is 
better matched to a non-orthogonal multiple-access backhaul: 
each receive sensor amplifies and forwards the received signal 
to the fusion center so that the signals transmitted by the 
receive sensors are superimposed at the fusion center (see, 
e.g., m, IH). 

Our specihc contributions are as follows: 

• CF: The joint optimization of the waveform and the quantiza¬ 
tion strategy is investigated for CF, with a focus on orthogonal- 
access backhaul. To reflect practical constraints, only stochas¬ 
tic channel state information (CSI) is assumed on the channel 
gains between target or clutter and the receive sensors. For 
an optimization objective, we adopt the information-theoretic 
criterion of the Bhattacharyya distance in order to account for 
the detection performance ISjl-lTl. 

• AF: The joint optimization of the waveform and the 
amplifying gains of the receive sensors is studied for AF, by 
concentrating on non-orthogonal multiple-access backhaul. We 
adopt the performance criterion and main assumptions of CF. 
Furthermore, we consider both instantaneous and stochastic 
CSI on the receive sensors-to-fusion center channels. 

Throughout, we assume tractable and well accepted models 
in order to gain insight into the problem at hand. With this 
insight gained, subsequent work may explore more detailed 
configurations. 


The rest of the paper is organized as follows. In Section 
im we present the signal model and cover the two types of 
backhaul links, namely orthogonal-access and non-orthogonal 
multiple-access backhaul. In Section [Hll after describing the 
CF backhaul transmission strategies and reviewing the optimal 
detectors, we present the optimization of the multistatic cloud 
radar system with CF. In Section IIVI we focus on the AF 
backhaul transmission, and optimize the system with both 
instantaneous and stochastic CSI under AF. Numerical results 
are provided in Section |V] and. Anally, conclusions are drawn 
in Section |VT] 

II. System Model 

Consider a multistatic cloud radar system consisting of a 
transmit element, N receive sensors, and a fusion center, or 
cloud processor, as illustrated in Fig. [T] The receive sensors 
communicate with the fusion center over an orthogonal-access 
backhaul or a non-orthogonal multiple-access backhaul. All 
the nodes are equipped with a single antenna, and the set of 
receive sensors is denoted Af = {l,...,iV}. 

The system aims to detect the presence of a single stationary 
target in a clutter held. To this end, each sensor receives a 
noisy version of the signal transmitted by the transmit element 
and reflected from the surveillance area, which is conveyed to 
the fusion center on the backhaul channels after either quan¬ 
tizing or amplifying the received signals as discussed below. It 
is assumed that perfect timing information is available at the 
fusion center, such that samples of the received signal may be 
associated with specific locations in some coordinate system. 
For such a location, and based on all the signals forwarded 
from the different receive sensors, the fusion center makes a 
decision about the presence of the target (see, e.g. Q-iini). 
Note that, as argued in a, the assumption of stationary target 
and scatterers can be regarded as a worst-case scenario for 
more general set-ups with non-zero Doppler. 

We consider a pulse compression radar in which the trans¬ 
mitted signal given in baseband form is 

K 

s{t) = ^ - (k- l)Tc), (1) 

k^l 

where (p{t) is, for example, a square root Nyquist with chip 
rate 1/Tc, so that {k — are orthonormal; and 

{^k}k=i ^ sequence of (deterministic) complex coefficients 

that modulate the waveform. The vector x = [xi ■ ■ ■ xk]"^ is 
referred to as waveform or code vector, on which we impose 
the transmit power constraint x^x < Pt- The design of the 
waveform x determines both target and clutter response, and 
thus has a key role in the performance of the radar system. 

The baseband signal received at the sensor n G Af, which 
is backscattered by a stationary target, can be expressed as 

rn{t) =hnS{t-Tn)+Cn{t)+Wn{t), (2) 

where is the random complex amplitude of the target 
return, which includes the effects of the channel and follows a 
Swerling I target-type model having a Rayleigh envelope, i.e., 
hn ^ CA/'(0, tT( „); Cn{t) represents the clutter component; 
Wn{t) is a Gaussian random process representing the signal- 
independent interference, which aggregates the contributions 
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of thermal noise, interference and jamming and is assumed 
to be correlated over time, as detailed below; and is the 
propagation delay for the path from the transmit element to 
the target and thereafter to the sensor n, which is assumed 
to satisfy the condition > KTc in order for the target 
to be detectable. The clutter component c„(f) consists of 
signal echoes generated by stationary point scatterers, whose 
echoes have independent return amplitudes and arrival times. 
Accordingly, the clutter component c„(f) is expressed as 

Nc 

^ ^ '^rL^v)t (3) 

v—1 

where is the number of point scatters; gn^v is the amplitude 
of the return from scatterer u; and „ is the propagation delay 
for the path from the transmit element to the scatterer v and 
to the sensor n, which satisfies the condition „ < KTc- 

After matched filtering of the received signal ^ with the 
impulse response </>*(—i), and after range-gating by sampling 
the output of the matched filter at the chip rate, the discrete¬ 
time signal at receive sensor n for n G Af can be written 
as 

'^n,k ~ “t” 9n^k “t“ '^n,k^ (4) 

where r„ ^ is the output of the matched filter at the receive 
sensor n sampled at time t = {k — l)Tc + Tn', the term gn = 

J2v=l9n,v'^{Tri - Tn,v) with T'(<) = <f>(T - {T)dT 

being the auto-correlation function of represents the con¬ 
tribution of clutter scatterers, which can be modeled, invoking 
the central limit theorem, as a zero mean Gaussian random 
variable with a given variance (see ||7] Appendix A]); 
and Wn^k is the fcth sample of Wn{t) after matched filtering at 
the sensor n. 

In vector notation, we can write (|4|i as 

fn — “f Cn (5) 

where we defined r„ = [r„ i • • • rn,K]'^, Sn — hnX and 
Cn — gnX', and the noise vector Wn — [wn,i ■■■ Wn,KY' 
follows a zero-mean Gaussian distribution with temporal cor¬ 
relation Vlyj^n, i-e., Wn CA/'(0,The variables /i„, gn 
and Wn for all n G Af, are assumed to be independent for 
different values of n under the assumption that the receive 
sensors are sufficiently separated ||6l. Moreover, their second- 
order statistics cr^ „ and are assumed to be known to 
the fusion center, for all n G Af, e.g., from prior measurements 
or prior information M, M- 

To summarize, the signal received at sensor n can be written 
as 

Ho ■rn=Cn+Wn, (6a) 

'Hi:rn=Sn+Cn+Wn, u G Af , (6b) 

where "Ho and 'Hi represent the hypotheses under which the 
target is absent or present, respectively. 

In the rest of this section, we detail the assumed model 
for both orthogonal-access and non-orthogonal multiple-access 
backhaul. 

Orthogonal-access Backhaul: For the orthogonal-access 
backhaul case, each receive sensor n is connected to the fusion 
center via an orthogonal link of limited capacity Cn bits per 
received sample. The capacity Cn is assumed to be known 


to the fusion center for 2 \\ n G Af and to change sufficiently 
slowly so as to enable the adaptation of the waveform and of 
the transmission strategy of the sensors to the values of the 
capacities C„ for all n G Af. 

Non-orthogonal Multiple-access Backhaul: For the non- 
orthogonal multiple-access backhaul, the signal received at 
the fusion center is the superposition of the signals sent by 
all receive sensors, where channels are subject to mutual 
interference. Accordingly, the received signal at the fusion 
center f = [fi • • • tkY' is given by 

N 

f = '^fntn+Z, ( 7 ) 

n—1 

where tn = [tn,i ■ ■ ■ is the signal sent by the receive 

sensor n on the backhaul to the fusion center; /„ is the 
complex-valued channel gain between the receive sensor n and 
the fusion center; and z — [zi ■ ■ ■ zk]'^ is the noise vector 
having a zero-mean Gaussian distribution with correlation 
matrix i.e., 2 ^ CAf Based on prior information 
or measurements, the second-order statistics of the channel 
gains between the target and the receive sensors, and of the 
noise terms, namely cr^„, a'^n^ ^w,n and are assumed to 
be known to the fusion center for all n G Af. The channel 
between receive sensors and fusion center / = [/i • • • /at]^ 
are also assumed to be known at the fusion center, via training 
and channel estimation. 

III. CF Backhaul Transmission 

In this section, we consider orthogonal-access backhaul and 
CF transmission. With CF, each receive sensor quantizes the 
received vector in (|6|l, and sends a quantized version of 
to the fusion center. Note that, since the receive sensor does 
not know whether the target is present or not, the quantizer 
cannot depend on the correct hypothesis Ho or Hi. In order to 
facilitate analysis and design, we follow the standard random 
coding approach of rate-distortion theory of modeling the 
effect of quantization by means of an additive quantization 
noise (see, e.g., El, El) as in 

rn=rn+qn, (8) 

where = [f„_i • • • is the quantized signal vector 

of r„; and ^ CAf {0,flq^n) is the quantization error vector, 
which is characterized by a covariance matrix Based 

on random coding arguments, while ([8]l holds on an average 
over randomly generated quantization codebooks, the results 
derived in this paper can be obtained by means of some (de¬ 
terministic) high-dimensional vector quantizer (see, e.g., El). 
For instance, as discussed in ll20l . a Gaussian quantization 
noise with any covariance „ can be realized in practice 
via a linear transform, obtained from the eigenvectors of 
followed by a multi-dimensional dithered lattice quantizer such 
as Trellis Coded Quantization (TCQ) 1211 . 

Based on (|8]l, the signal received at the fusion center from 
receive sensor n is given as 

Ho . f'n — Cn “t" Wn Qm 
Hi : rn=Sn+Cn+Wn+ Qn- 
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As further elaborated in the following, the covariance matrix 
Qg jT, determines the bit rate required for backhaul communi¬ 
cation between the receive sensor n and the fusion center iflTl . 
m and is subject to design. 

To set the model (|9|l in a more convenient form, the signal 
received at the fusion center is whitened with respect to the 
overall additive noise c„ + Wn + Qn, and the returns from all 
sensors are collected, leading to the model 

■Hi : j/-CAr(0,T>5£> + /), ^ 

whem y=[yj' ■■■ 2/^]^, J/n = D^fn, is the whitening 
matrix associated with the receive sensor n and is given by 
Dn = D is the block diagonal 

matrix D = diag{Z)i,Djv}, and S is the block diago¬ 
nal matrix S = The detection 

problem formulated in ( fTUl l has the standard Neyman-Pearson 
solution given by the test 

l~i\ 

y^Ty%u, (11) 

Ho 

where we have defined T = DSD {DSD +I)-\ and the 
threshold i/ is set based on the tolerated false alarm probability 

m- 

In the rest of this section, we aim to hnd the optimum 
code vector x and quantization error covariance matrices fig „ 
in for given backhaul capacity constraints Cn, for all 
n G Af. Before we proceed, for reference, we first discuss 
the standard distributed detection approach that combines hard 
local decisions at the receive sensors and a majority-rule 
detection at the fusion center (see, e.g., US, El). 

A. Distributed Detection 

Here, we describe the standard distributed detection ap¬ 
proach applied to multistatic radar system (see, e.g., E3l . 
El)- With this approach, each receive sensor n makes its own 

■Hi 

decision based on the likelihood test given by y^Tnyn I 7n, 

-Ho 

where 7 „ is the threshold for receive sensor n, which is 
calculated based on the tolerated false alarm probability E2l . 
and we have dehned T„ = DnSnDn{DnSnDn + I)~^ 
with j/„ = -D„r„, Dn = {aln^x^ + and 

Sn — for all n G Af. The receive sensors transmit 

the obtained one-bit hard decision to the fusion center. Note 
that this scheme is feasible as long as the backhaul capacity 
available for each receive sensor-to-fusion center channel is 
larger than or equal to 1/H bits/sample, i.e., Cn > 1/H, for 
n G Af. The fusion center decides on the target’s presence 
based on the majority rule: if the number of receive sensors 
k that decide for Ho satisfies k > N/2, the fusion center 
chooses Ho, and vice versa if fc < N/2. 


B. Performance Metrics and Constraints 

To start the analysis of the cloud radar system, we discuss 
the criterion that is adopted to account for the detection perfor¬ 
mance, namely the Bhattacharyya distance and the approach 


used to model the effect of the quantizers at the receive 
sensors. 

Bhattacharyya Distance: For two zero-mean Gaussian 
distributions with covariance matrix of Ei and E 2 , the Bhat¬ 
tacharyya distance B is given by ||5l 


H = log 


( |0.5(Si-fS2)| ^ 

vIsIW )' 


( 12 ) 


Therefore, for the signal model ®, the Bhattacharyya distance 
between the distributions under the two hypotheses can be 
calculated as 


B{x,Slg) = log 


\I + 0.5DSD\ 


v/|/ + DSD 

N 

— 'y ^. Bn{x^Clqn) 

1 + O.bX 


n—1 
N 


= ^log 


v/l + An 


(13) 


where we have made explicit the dependence on x and fig 
fig collects all the covariance matrices of quantization noise 
and is given as fig = {flg,n}neH; we have defined 

\n — A' flu),n A ^q^n^ X. (14) 

We observe that (fOl l is valid under the assumption that the 
effect of the quantizers can be well approximated by additive 
Gaussian noise as per ®. This is discussed next. 

Quantization: From rate-distortion theory, a vector quan¬ 
tizer exists that is able to realize the additive quantization noise 
model ([8]), when operating over a sufficiently large number 
of measurement vectors as long as the capacity C„ is 
no smaller than the mumal information /(r„;f„)/H lfT9l . 
For example, a dithered lattice vector quantizer achieves this 
result EOl . These considerations motivate the selection of the 
mutual information as a measure of the backhaul 

rate required for the transmission to the fusion center. 

While the mutual information depends on the 

actual hypothesis Ho or Hi, it is easy to see that /(r„;r„) 
is larger under hypothesis Hi. Based on this, the mutual 
information evaluated under Hi is adopted here as 

the measure of the bit rate required between receive sensor 
n and the fusion center. This can be easily calculated as 
I{rn\ Tn) = In{x, Clq^n) by usiug the expression of the mutual 
information for multivariate Gaussian distribution (see, e.g., 
ifT^ ) with 

(X, fig^Tl) = log |/ -f (fig^Tl) 

-flog (1 -f {a/n A aln)x^{^w,nA^q,n)~^x) , (15) 


where again we have made explicit the dependence of mutual 
information on x and fig,ra¬ 
in the following, we formulate and solve the problem of 
jointly optimizing the Bhattacharyya distance criterion over the 
waveform x at the transmit element and over the covariance 
matrices fig of the quantizers at the receive sensors in Section 
IIII-CI and in Section IIII-Dl respectively. 
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C. Problem Formulation 

The problem of maximizing the Bhattacharyya distance in 
(foi l over the waveform x and the covariance matrices Clq 
under the backhaul capacity constraints is stated as 

N 

minimize B{x,Vtq) = (16a) 

n—1 

S.t. In{x,ilq^n) < KCn = Cn, UG Af, (16b) 

x^x < Pt, (16c) 

^q,n ^0, n G Af, (16d) 

where we have formulated the problem as the minimization 
of the negative distance B{x,Slq) = J2n=i ^n(x,^q,n), with 
B{x,Qq) = -B{x,Slq) and 

following the standard convention in The power of the 
waveform x is constrained not to exceed a prescribed value 
of transmit power Pt- We observe that the constraint (I16bl i 
ensures that the transmission rate with K chips between each 
receive sensor and the fusion center is smaller than Cn, 
according to the adopted information-theoretic metric. Note 
also that the problem ( fThl l is not a convex program, since 
the objective function (I16ab and the constraints (I16bb are not 
convex. 

D. Proposed Algorithm 

Since both functions B„(x,Qq^) and Xn(x,iiq„) in (fThb are 
non-convex in x and the optimization problem (fThb is 

not convex, and hence it is difficult to solve. To obtain a locally 
optimal solution, we approach the joint optimization of x and 
ilq in (fThl l via successive convex approximations. Specifically, 
in an outer loop. Block Coordinate Descent (BCD) is applied 
to update x and Qq one at a time, while an inner loop 
implemented via Majorization-Minimization (MM) solves the 
optimization of x and flq separately. This approach was first 
introduced in ll^ for a sum-capacity backhaul constraint. By 
the properties of MM (see, e.g., Il27l . Il28l l. the algorithm pro¬ 
vides a sequence of feasible solutions with non-increasing cost 
function, which guarantees convergence of the cost function. 
Note that, due to the non-convexity of the problem, no claim 
of convergence to a local or global optimum is made here. 

At the ith iteration of the outer loop, the optimum waveform 
x6) is obtained by solving (fThb for matrices Clq = 
obtained at the previous iteration; subsequently, the matrices 
are calculated by solving (fThb with x = These two 
separate optimizations are carried out by the MM method, 
which, as described in Appendix A, requires the solution of 
a quadratically constrained quadratic programs (QCQP). The 
proposed algorithm coupling BCD and MM to solve problem 
(fThb . is summarized in Table Algorithm [T] In Algorithm [T] we 
use the superscript i to identify the iterations of the outer loop, 
and the superscript j as the index of the inner iteration of the 
MM method (e.g., indicates the waveform optimized at 
the jth iteration of the inner loop of the MM method and the 
ith iteration of the outer loop). In Appendix A, we present the 
MM steps and the overall proposed algorithm in detail. 

The complexity of Algorithm [T| by using standard convex 
optimization tools is polynomial in K and N since, at each 


outer iteration, MM requires to solve the problems (Ia3T i 
and whose sizes of the optimization domains are K 

and NK"^, and numbers of constraints are iV -f 1 and 27V, 
respectively ll25l. 1291. 


Algorithm 1 Joint optimization of waveform and quantization 
noise covariances (fThb 

Initialization (outer loop); Initialize x^°^ G ^ 

0 and set 7 = 0. 

Repeat (BCD method) 

% \ 

Initialization (inner loop): Initialize = 

x(®“i) and set j = 0. 

Repeat (MM method for x^*)) 

j ^ j + 1 

Find x^'^’P by solving the problem (IA.3b with 

Until a convergence criterion is satisfied. 

Update x^*) ^ x^^’P 

Initialization (inner loop); Initialize = 

and set j = 0. 

Repeat (MM method for 

^ ^ t \ 

Find by solving the problem ( IA.6b with 

X = x^^. 

Until a convergence criterion is satisfied. 

Update Clq'^ g- Vlq'^^ 

Until a convergence criterion is satisfied. 

Solution: x •(— x'P and Q,q g- 


IV. AF Backhaul Transmission 

In this section, we consider AF transmission on a non- 
orthogonal multiple-access backhaul. With AF, sensor n G Af 
amplifies the received signal r„ in (|6ll and then forwards the 
amplified signal tn = a.ni'n to the fusion center, where an is 
is the amplification coefficient at the receive sensor n. From 
(01, the fusion center is faced with the following detection 
hypothesis problem 

N N 

1~Lq . T — ^ ^ fnin ^ ^ ^ fnC^n (^n “t" ^n) 

n—1 n—1 

N N 

Ph: r = fntn jnOin (s„ + Cn + Wn) + Z. 

n—1 n—1 

(17) 

The variables /i„, pn, Wn, fn and z, for all n G Af, are 
assumed to be mutually independent. Since only the second- 
order statistics of the channel gains hn, n G Af, are known 
to the receive sensors and the fusion center, no coherent 
gains may be achieved by optimizing the amplifying gains, 
and hence one can focus, without loss of optimality, only 
on the receive sensors’ power gains p = [pi • • -pnY, with 
Pn = \an\^, for n gN. 

As in the CF backhaul transmission in Section |III] we can 
write the hypotheses (fTTl i in a standard form by whitening 
the signal received at the fusion center, and consequently 
the detection problem can be expressed as (O, where we 
have redefined y = DP, D = {Yjn=ii\fn?Priori,n^x^ + 
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\fn\‘^Pn^w,n) + IS the whitening hlter with respect 

to the overall additive noise fnCtnicn + Wn) + -z; and 

S = '^^e Correlation matrix of the 

desired signal part. Accordingly, the detection problem has 
the standard estimator-correlator solution given by the test 
in (fTTl i. In the rest of this section, we seek to optimize the 
detection performance with respect to the waveform x and 
the power gains p, under power constraints on the transmit 
element and receive sensors. As done above, we adopt the 
Bhattacharyya distance as the performance metric. As per 
(fT2l l. the Bhattacharyya distance between the distributions (fTTl i 
of the signals received at the fusion center under the two 
hypotheses T-Lq and T-Li can be calculated as 




log 


log 


\I + 0.5DSD\\ 
^\I + DSD\ ) 


1 + 0.5A\ 

7r+Tj’ 


(18) 


where A = + (/ 0 Ik)^ (P 0 

lK)^w{f ^ Ik) + St = ,crt^_tv} 

and Sc = diag{CTc 1 ,..., CTc wi *-1*6 diagonal matrices 
whose components are the second-order statistics of channel 
amplitudes of target return and clutter, respectively; — 
diag{f 2 u, 4 ,... jiltu tv} S ^ block diagonal ma¬ 

trix containing all the noise covariance matrices at the receive 
sensors; and P = diag{p} S ji^xN (jjagonal matrix that 
contains the receive sensors’ power gains. Note that we have 
made explicit the dependence of the Bhattacharyya distance 
B{x,p;f) on the channels / at the fusion center, as well as 
on the waveform x and the receive sensors’ power gains p. 


A. Short-Term Adaptive Design 

We first consider the case in which design of the wave¬ 
form X and of the receive sensors’ gains p depends on 
the instantaneous gain of the CSI of the receive sensors-to- 
fusion center channels /. Note that this design requires to 
modify the solution vector {x,p) at the time scale at which 
the channel vector / varies, hence entailing a potentially 
large feedback overhead from the fusion center to the receive 
sensors and the transmit element. The problem of maximizing 
the Bhattacharyya distance (fTsT i over the waveform x and the 
power gains p under the power constraints for transmit element 
and receive sensors, is stated as 


minimize 

x,p 

B{x,p-,f) 

(19a) 

s.t. 

x^x < Pt, 

(19b) 


< Pr, 

(19c) 


Pn>0, n G Af, 

(19d) 


where we have defined B{x,p;f) = —B{x,p;f) to formulate 
the problem as the minimization of the negative Bhattacharyya 
distance B{x,p;f). We observe that the problem ( fT9b may be 
easily modihed to include individual power constraints at the 
receive sensors, but this is not further explored here. Moreover, 
the problem (fT^ is not a convex program, since the objective 
function ( I19al i is not convex. 


We propose an algorithm to solve the optimization problem 
(ESll. As in Section IIII-DI due to the difficulty of obtaining 
a global optimal solution, we develop a descent algorithm, 
and adopt the BCD method coupled with MM. The proposed 
algorithm is summarized in Table Algorithm |2] and further 
detailed in Appendix B. The complexity of the Algorithmic 
by using standard convex optimization tool is polynomial in 
K and N since, at each outer iteration, MM requires to solve 
the problems (IB.21 l and (IB.41 l. whose sizes of the optimization 
domains are K and N, and numbers of constraints are 1 and 
N 1, respectively Il25l . 1291 . 


Algorithm 2 Short-term adaptive design of waveform and 
ampliher gain (fT^ 

Initialization (outer loop); Initialize x^^'> € p^^'> ^ 

0 and set i = 0. 

Repeat (BCD method) 

i i — i T 1 

Initialization (inner loop); Initialize = 

and set j = 0. 

Repeat (MM method for x^*^) 

J ^ J + 1 

Find by solving the problem ( IB.21 i with 

p — p^'‘~^\ 

Until a convergence criterion is satisfied. 

Update x^*) ^ 

Initialization (inner loop); Initialize = 

p(*-i) and set j = 0. 

Repeat (MM method for p^*)) 

j ^ j + 1 

Find by solving the problem ( IB.41 i with 

X = x^^\ 

Until a convergence criterion is satisfied. 

Update ^ 

Until a convergence criterion is satisfied. 

Solution; x ^ x^*) and p ^ p(*) 


B. Long-Term Adaptive Design 

Here, in order to avoid the possibly excessive feedback 
overhead between fusion center and the transmit element 
and receive sensors of the short-term adaptive solution, we 
adopt the average Bhattacharyya distance, as the performance 
criterion, where the average is taken with respect to the 
distribution of the receive sensors-to-fusion center channels 
/. In this way, the waveform x and receive sensors’ gains 
p have to be updated only at the time scale at which the 
statistics of channels and noise terms vary. Then, the problem 
for the long-term adaptive design is formulated from problem 
(fT^ by substituting the objective function B(x,p;f) with 
F;/[e(x,p;/)], yielding 

minimize Ef [S(x,p;/)] (20a) 

s.t. ([Hbll - (dldll. (20b) 

Note that the problem (l20l i is a stochastic program with a 
non-convex objective function (I20ab . 

Since the stochastic program ( l20b has a non-convex ob¬ 
jective function, we apply the stochastic successive upper- 
bound minimization method (SSUM) fSO), which minimizes at 
each step an approximate ensemble average of a locally tight 
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Fig. 2. Bhattacharyya distance versus the backhaul capacity Cn = C, n £ J\f 
for CF backhaul transmission, with Pt = 10 dB, K = 13, = 3, cr^ ^ = 1, 

0-2 1 = 0.125, 0-2 2 = 0.25, 0-2 3 = 0.5 and = (1 - 0.12n)l»--’'l 

for n £ M. 

upper bound of the cost function. Specifically, we develop a 
BCD scheme similar to the one detailed in Table Algorithm 
|2] that uses SSUM in lieu of the MM scheme. Details are 
provided in Appendix C. The final algorithm for long-term 
adaptive design can be summarized as in Table Algorithm[2]by 
substituting (IB.2b and (IB.4b with (1C.lb and (1C.2b . respectively. 
Convergence of the SSUM algorithm is proved in If^ and 
the algorithm guarantees feasible iterates. The complexity of 
the proposed algorithm by using standard convex optimization 
tool is polynomial in K and N since, at each outer iteration, 
SSUM requires to solve the problems (1C. lb and (1C.2b . whose 
sizes of the optimization domains are K and N, and numbers 
of constraints are 1 and N + 1, respectively ll25l . Il29l . 

V. Numerical Results 

In the following, the performance of the proposed algo¬ 
rithms that perform joint optimization of the waveform x and 
of the quantization noise covariance matrices Clq for the CF, 
and of the waveform x and of the power gains p for AF, 
are investigated via numerical results in Section IV-AI and in 
Section IV-BI respectively. Throughout, we set the length of 
the waveform to AT = 13 and the variances of the target 
amplitudes as = 1 for n G Af. For reference, we 
consider a baseline waveform with Barker code of length 13, 
i.e., 6 i 3 = [1 1 1 1 1 - 1 -111 - 11-1 1 ]^. 
Moreover, unless stated otherwise, we model the noise with 
covariance matrices [Clw,n]i,j = (1 — 0 . 12 n)l®“'^l and = 

(1 - 0 . 6 )I*-Jl as in Q, hence accounting for temporally 
correlated interference. The channel coefficients /„ have unit 
variance, i.e., = 1 . 

A. CF Backhaul Transmission 

In this section, the performance of the proposed joint 
optimization of the waveform x and of the quantization noise 
covariance matrices in Section Hill is verified via numerical 
results. Note that some limited results for a sum-backhaul 
constraint were presented in 12^ . For reference, we consider 
the performance of the upper bound obtained with infinite 



Fig. 3. Probability of detection Pd versus the backhaul capacity Cn = C 
for CF backhaul transmission, n £ J\f, with Pt = 10 dB, K = 13, N = 3, 
= 1, = 0.125, = 0-25, = 0.5, [nn,,n]i,j = (1 - 



Fig. 4. ROC curves for CF backhaul transmission with Pt = 10 dB, 
K = 13, N = 3, = 1, 0-23 = 0.125, = 0.25, 0-23 = o.5, 

[ilw,n]i,j = (1 — 0.12n)b“.tl and Cn = (7 = 5 for n € Af. 

capacity backhaul links, distributed detection using the Barker 
waveform (see. Section IIII-Ab . and the following strategies; 
(0 No optimization (No opt): Set x = Pt/K bi^ and 
^q,n = cl, for n G Af, where e is a constant that is found 
by satisfying the constraint (I16bb with equality; (ii) Waveform 
optimization (Waveform opt.) : Optimize the waveform x by 
using the algorithm in JT), which is given in Algorithm [T] 
by setting = 0 for n G Af, and set = el, for 

n G Af, as explained above; (Hi) Quantization noise optimiza¬ 
tion (Quantization opt): Optimize the covariance matrices ilq 
as per Algorithm [T] with x — QPr/Kbi^. In the following, 
we set the number of receive sensors, the transmit power and 
the variance of the clutter amplitudes as V = 3, Pt = 10 dB, 
(Tg 1 = 0.125, (Tg 2 = 0.25 and cr^ g = 0.5, respectively. Also, 
the backhaul rate constraints Cn are assumed to be equal, i.e., 
Cn = C for all n G Af. 

In Fig. |2] the Bhattacharyya distance is plotted versus the 
available backhaul capacity C. For intermediate and large 
values of C, the proposed joint optimization of waveform 
and quantization noise is seen to be significantly beneficial 
over all separate optimization strategies. In order to study the 
actual detection performance and validate the results in Fig. 
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(c) Optimal waveform and quantization noise with low-frequency (d) Optimal waveform and quantization noise with high-frequency 
interference interference 


Fig. 5. Comparison of the energy/power spectral densities of the waveforms obtained with a Barker code (Barker waveform) and with an optimal code x 
(Optimal waveform), and of optimal quantization noise obtained by Algorithm ^ when Pt = 10 dB, K = 13, TV = 3, = 1, ^ = 0.125, 

2 = 0.25, cTc 3 “ = C = 5 for n C (a) and (c) consider receive sensors with low-frequency interference having temporal con'elation 

[^w,n]i,j = (1 — 0.12n)l*“-?l, (b) and (d) consider receive sensors with high-frequency interference having temporal correlation [Qw,n]i,j = (—1 + 

0.12n)l*-jl. 


[2l Fig- E] shows the detection probability Pd as a function 
of the available backhaul capacity C when the false alarm 
probability is Pfa = 0.01. The curve was evaluated via Monte 
Carlo simulations by implementing the optimum test detector 
(dUi. We also implemented the distributed detection scheme 
described in Section IIII-AI by setting the threshold 7 „ to be 
equal for n G Af for simplicity. It can be noted that the relative 
gains predicted by the Bhattacharyya distance criterion in Fig. 
|2] are consistent with the performance shown in Fig. 0 More¬ 
over, for small values of C, distributed detection outperforms 
cloud detection due to the performance degradation caused 
by the large quantization noise on the cloud-based schemes. 
However, as the available backhaul capacity C increases, the 
cloud detection approach considerably outperforms distributed 
detection. 

Fig. |4] plots the Receiving Operating Characteristic (ROC), 
i.e., the detection probability versus false alarm probability 
Pfa, for (7 = 5. It is conhrmed that the proposed joint opti¬ 
mization method provides remarkable gains over all separate 
optimization schemes as well as over the distributed detection 


approach. For instance, for Pfa = 0.01, joint optimization 
yields Pd = 0.7251, while waveform optimization only yields 
Pd = 0.4556. 

Fig. |5] shows the energy/power spectral density functions 
of the waveform with Barker code (Barker waveform) and 
with optimal code x (Optimal waveform), and of optimal 
quantization noise {qn}n=i obtained by Algorithm [T] when 
a square root Nyquist chip waveform (j>(t) with duration Tc 
is adopted, and Cn = C = 5 for n G Af. We consider two 
types of interference at the receive sensors, namely (a) low- 
frequency interference with temporal correlation [Ctaj^n]i,j = 
(1 — 0.12n)l*“'^l which has a single spectral peak at zero 
frequency; and (b) high-frequency interference with temporal 
correlation = (—1-f 0.1271)1*“-^I, having a minimum 

at zero frequency. It is observed in Fig. |5(c)| and Fig. |5(d)| 
that the spectrum of the optimal waveform concentrates the 
transmitted energy at frequencies for which the interference 
power is less pronounced, while the spectrum of the quantiza¬ 
tion noise concentrates at frequencies and sensors for which 
the interference power is more pronounced. 
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Fig. 6 . Bhattacharyya distance versus the transmit element’s power Pt for 
AF backhaul transmission with Pn = 10 dB, K = 13 , N = 3 , = 

1, al„ = 1, 0-2 ^ = 0.25, = 0.5, cr^g = 1 , = (1 - 

0.12n)0“fl and ['nz]ij' = (1 — 0.6)1*“^! for n € Af. 

B. AF Backhaul Transmission 

In this section, we evaluate the performance of the proposed 
algorithms that perform the joint optimization of the waveform 
X and of the amplifying power gains p for the short-term 
(Section FlV-Ab and long-term (Section FlV-Bb adaptive designs. 
For reference, we consider the following schemes; (i) No opt.: 
Set X = ^Pt/K bi 3 and p = Pr/NIn; (ii) Waveform opt: 
Optimize the waveform x as per Algorithm (with (1C. Il l in 
lieu of (IR2 T i for the long-term adaptive design) with p — 
Pr/N1n', and (Hi) Gain optimization (Gain opt): Optimize 
the gains p as per Algorithm [2] (with (1C.21) in lieu of (IB .41 ) for 
the long-term adaptive design) with x = Pt/K bi^. We set 
the total receive sensors’ power as Pr = 10 dB. Note that the 
upper bound with ideal backhaul is far from the performance 
achieved with AF over a non-orthogonal backhaul even for 
large sensors’ power Pr, and it is hence not shown here. The 
gap between the AF performance and the upper bound is due to 
the fact that, in order to obtain an ideal backhaul, one needs to 
code across long block lengths whereas AF operates on block 
length of size equal to the waveform K (here K = 13). 

Fig. | 6 ] shows the Bhattacharyya distance as a function of 
the transmit element’s power Pt, with N = 3, ^ = 0.25, 

(7^2 = 0-5 Pc 3 — 1- Fof small values of Pt, optimizing 
the waveform is more advantageous than optimizing the am¬ 
plifying gains, due to the fact that performance is limited by 
the transmit element-to-receive sensors connection. In contrast, 
for intermediate and large values of Pt, the optimization of 
the receive sensors’ gains is to be preferred, since the perfor¬ 
mance becomes limited by the channels between the receive 
sensors and the fusion center. Joint optimization significantly 
outperforms all other schemes, except in the very low- and 
large-power regimes, in which, as discussed, the performance 
is limited by either the transmit element-to-receive sensors 
or the receive sensors-to-fusion center channels. In addition, 
we observe that the long-term adaptive scheme loses about 
30% in terms of the Bhattacharyya distance with respect to 
the short-term adaptive design in the high SNR regime. The 
results in Fig. | 6 ] can be interpreted by noting that the joint 



Fig. 7. Bhattacharyya distance versus the number receive sensors N for AF 
backhaul transmission with Pt = 5 dB, Pr = 10 dB, K = 13, = 1, 

= T <72 g = 1, <72 2 = 0.9, 723 = 0.75, = 0.5, <72 ^ = 0.35, 

0-2 g = 0.25 and 0-2 ^ = 0.125, o-2 g = 0.05, = (l-0.12n)l‘-^l 

and = (1 — 0.6)h“7l for n S A/”. 



Fig. 8 . ROC curves for AF backhaul transmission with Pt = 5 dB, Pr = 
10 dB, A' = 13, Af = 3, O-^^ = 1, o-2„ = 1, 0-2^ = 0.25, = 0.5, 

0-23 = 1 , = (1 - 0 . 12 n)h-ll and = (1 — 0 . 6 )h-jl for 

nkN. 

optimization seeks to design the transmitted signal x such that 
it reduces the power transmitted at the frequencies in which the 
receive sensors observe the largest interference, while, at the 
same time, allocating more power to receive sensors suffering 
from less interference and, with the short-term adaptive design, 
having better channels to the fusion center. 

In Fig. |7] the Bhattacharyya distance is plotted versus the 
number receive sensors N with Pt = 5 dB, = 1, 
<2 = 0.9, <3 = 0.75, <4 = 0.5, <5 = 0.35, <6 = 0.25, 

7 = 0.125 and = 0.05. Optimizing the receive sensors’ 
power gains is seen to be especially benehcial at large N, due 
to the ability to allocate more power to the receive sensors 
in better condition in terms of interference and channels 
to the fusion center. For instance, even with the long-term 
adaptive design, optimizing the receive sensors’ power gains 
outperforms waveform optimization with short-term adaptive 
design for sufficiently large N. 

Fig. [ 8 ] plots the ROC curves with Pt = 5 dB, N = 3, 
cr'^i = 0.25, 2 = 0.5 and cr^ 3 = 1. The curve was evaluated 
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via Monte Carlo simulations by implementing the optimum 
test detector (fTTT i as discussed in Section|IV] It can be observed 
that the gains observed in the previous figures directly translate 
into a better ROC performance of joint optimization. Note also 
that power gain optimization is seen to be advantageous due 
to sufficient value of Pt as predicted based on Fig. |6] 

VI. Concluding Remarks 

We have studied a multistatic cloud radar system, where 
the receive sensors and fusion center are connected via an 
orthogonal-access backhaul or a non-orthogonal multiple- 
access backhaul channel. In the former case, each receive 
sensor quantizes and forwards the signal sent by transmit 
element to a fusion center following a compress-and-forward 
protocol, while amplify-and-forward of the received signal 
is carried out over the multiple-access backhaul. The fusion 
center collects the signals from all the receive sensors and 
determines the target’s presence or absence. We have in¬ 
vestigated the joint optimization of waveform and backhaul 
transmission so as to maximize the detection performance. 
As the performance metric, we adopted the Bhattacharyya 
distance and the proposed algorithmic solutions were based on 
successive convex approximations. Overall, joint optimization 
was seen to have remarkable gains over the standard separate 
optimization of waveform and backhaul transmission. More¬ 
over, cloud processing is found to outperform the standard 
distributed detection approach as long as the backhaul capacity 
is large enough. 

Appendix A; Details oe CF Optimization 
A. Review of MM Method 

We start by reviewing the MM method. For a non-convex 
function f{f) of a generic variable t, which may appear 
either in the cost function or among the constraints, the MM 
method substitutes at the Zth iteration, a convex approximation 
of f(t), such that the global upper bound property 
> f{t) is satisfied for all t in the domain, along 
with the local tightness condition 

These properties guarantee the feasibility of all iterates and 
the descent property that the object function does not increase 
along the iterations. 


B. Details of the Proposed Algorithm\I\ 

In the following, we discuss the application of the MM 
method to perform optimizations over x and fig in Algorithm 
[1] respectively. 

Optimization over x: Here, the goal is to obtain the optimal 
value of for problem (fT^ given O, = To this end, 

we apply the MM method. Specifically, at the jth iteration of 
the MM method and the ith iteration of the outer loop, the 
MM method solves a QCQP and obtains a solution x*^*by 
substituting the non-convex objective function B{x,ilq) with 
a tight upper bound around the current 

iterate x^'‘’^~^\ This bound is obtained by linearizing the 
difference-of-convex functions in B{x,ilq) via the first-order 


Taylor approximation ||27| . which follows the same steps as 
in 10 eq. (34) and (50) in Section IV], and is given by 


N 




N 


= E' 

n—l 

where 

A(bf-i) = 




X {^w,n T ^q^n) ttC 


(A.l) 


Pn 


1 + PnV' 


0.57„ 


+ Pn{^ + 0.57„) 

Pn 


1 + a: 


(ij-l) 


1 + PnU 


b.i-i) 




2p (1 + 0.57„) 


1 + PnUn'^ (1 + 0 . 57 „) 

+2^„ (1 + 0.57„)) ; 

H 




-ixbj-1); 




1 + PnVn 

with Pn = of. n ^nd 7 „ = a^.^/ Pn- A bound with the desired 
property can also be easily derived for I„(x,Og „) by using 
the inequality log(l-l-f) < log(l-|-f^^^) + l/(l-l-f^*^)(f — 
for t = {al„ + of „)x^{ilyj^ri + ^q,n)~^x, leading to 

Uf;{x,nq^„\x^^’^~^'>) = log|7 -b 


+ log(l + {(^c,n + ^t,n) 


(A.2) 


At the jth iteration of the MM method and the zth outer loop, 
we evaluate the new iterate by solving the following 

QCQP problem 

x(®J) ^ argmin 7f®(x,(A.3a) 

X 

s.t. W^(x,^2,,„|x(*’^-l)) < C„, neV, (A.3b) 
x^x < Pt- (A.3c) 


The MM method obtains the solution x*^®^ for the fth iteration 
of the outer loop by solving the problem (IA.3I) iteratively over 
j until a convergence criterion is satisfied. 

Optimization over Clqi In this part, we consider the 
optimization of matrices for a given x = x*^®\ Sim¬ 
ilar to the optimization over x*^®\ we use upper bounds of 
B{x,flq) and I„(x,Og „) for optimization. First, by rewriting 
I„(x,Oq,„) as J„(x,Og,n) = log + {o^^ + af jxx^ + 
7iiu,n| — log|Og ,i|, we obtain difference-of-convex functions 
with respect to Clq^n- Then, by linearizing negative convex 
component via its first-order Taylor approximation, upper 
bounds and with the 

desired properties of MM method are derived for functions 
I„(x,f2g_„) and B{x,Ctq), respectively, as follows: 
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= log + {(jln + 

- log |0,,„| + tr I + {al^ + (jIJxx^ 

+n^,„)-' I (a.4) 

and 

N 

n—1 
N 

= ^ - log KO.ScTt^^^ + crln)xx^ + + f^g,n| 

n—1 

+0.5tr I + al^)xx^ + 
xfiq.n} + 0.5tl- I (al „^XX^ + ^w,n + 

(A.5) 

The jth iteration of the MM method then evaluates the ma¬ 
trices }n€j\r by solving the following convex 

optimization problem 


where 


0.57 P 

1 + AbO-1) _|_ ’ 

^ ( 2 ^( 1 + 0 . 57 ) 

\l + / 3 j/bo' 1) ( 1 -p 0.57) 

+ 2/3 (1 + 0.57)) {{f (g) Ik )^ {P (g) Ik ) 




/3 = /^PS,/; 

_ /^PSt/ 

^ ^ ' H 

((/ 0 Ij,f (p ® Ij,) 


A(*o-i) = 7 - 


7 

1 +/3y(*■•?-1) ■ 


At the jth iteration of the MM method and the ith outer loop, 
we evaluate the new iterate xl*’-!! by solving the following 
QCQP problem 

^ argmin U{x,p] f\x^'''^~^'^) (B.2a) 

X 


s.t. x^x < Pt- (B.2b) 


^ argminW®(x,fi,|f 2 lbJ-i)) (A. 6 a) 

fig 

S.t. < C„, n&M, (A. 6 b) 

^q,n ^0, n G Af. (A. 6 c) 

By repeating the procedure (IA. 6 I 1 over j until the convergence 
is attained, the solution Qq '^ is obtained for the ith outer loop. 


Appendix B: Details of AF Short-term Adaptive 
Design 

A. Optimization over x 

Here, the goal is to optimize the objective function ( fT^ over 
the waveformxl®! given the gainsp = For this purpose, 

we apply the MM method. Specifically, at the jth iteration of 
the MM method and the ith iteration of the outer loop, the MM 
method solves a convex QCQP and obtains a solution xl* -^! 
by substituting the non-convex objective function B{x,p;f) 
with a tight upper bound W(x,p; /|xl*A-i)) around the current 
iterate This bound is obtained by following the same 

steps as in Appendix lA-BI and is given by 


The MM method obtains the solution x*^®^ for the iih iteration 
of the outer loop by solving the problem (IB .2b iteratively over 
j until a convergence criterion is satisfied. 

B. Optimization over p 

We consider now the optimization of the gains when 
the waveform x = x^*^ is given. Similar to the optimization 
over in the previous section, we also use the MM method 
for the optimization over p. Towards this goal, we obtain 
the upper bound ///(x,p;of the objective function 
B{x,p]f) around the current iterate This bound is 

derived by linearizing the difference-of-convex functions via 
the first-order Taylor approximation lIZTll . The bound can then 
be obtained in ( IB.3b at the top of the next page. Then, the new 
iterate pb'i) at the jth iteration of the MM method and the 
/th iteration of the outer loop can be obtained by solving the 
following optimization problem: 

p(*j) ^ argmin U{x,p; f\p^'^’^~^'>) (B.4a) 

p 

s.t. I'^p < Pr, (B.4b) 

p„ > 0, n G Af. (B.4c) 

By repeating the procedure (IB .4b over j until a convergence 
criterion is satisfied, the solution p^*) is determined for the /th 
outer loop. 


Kix,p-f\x^^’^ 

= ((/ ® /k)" (P ® Ik ) (/ ® Ik ) 

(B.l) 


C. Summary of the Proposed Algorithm |2] 

In summary, in order to solve problem ( fT9b , we propose 
an algorithm (described in Table Algorithm [21l that alternates 
between the optimization over x, described in Appendix IB-AI 
and the optimization over p, discussed in Appendix IB-BI In 
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i^{x,p;f\p 




= - In 


/"P (0.5Et + E,) /xx" + (/ ® iKf (P ® /k) (/ ® /k) + 


+0.5tr I (/^p(*.^-i) (Si + Ee) /xx^ + (/ 0 iKf (p(*’^-i) 0 7 k) (/ 0 7k) + 0.) 

X (/^P (Ei + E,) fxx^ + (/ 0 7k)^ (P 0 7k) (/ 0 7k)) } 

+0.5tr I (/^p(*’^-i)E,/xx^ + (/ 0 iK f (p(*’^-i) 0 7k) (/ 0 7k) + J^z) 

X (/^PE,/xx^ + (/ 0 7k)^ (P 0 7k) (/ 0 7k)) } ■ 


(B.3) 


particular, at the ith iteration of the outer loop, the iterate 
x^*^ is obtained by solving a sequence of convex problems 
(Appendix IB-Ab via the MM method for a fixed p — 

Then, the iterate is found by solving a sequence of convex 
problems (Appendix IB-BI) via the MM method with x = x^*^ 
attained in the previous step. According to the the properties 
of the MM method IIZTl . Il28l . the proposed scheme yields 
feasible iterates and a non-increasing objective function along 
the outer and inner iterations, hence ensuring convergence of 
the cost function. 

Appendix C; Details oe AF Long-term Adaptive 
Design 

A. Optimization over x 

Following the SSUM scheme, at the jth inner iteration and 
the zth outer iteration, we optimize the waveform given 
p — p^^~^'> by solving the following convex problem 

1 ^ 

x(^J)^argmin - ^7((')(x,p;|x(*’'-i)) (C.la) 

® ^ 1=1 

s.t. x^x < Pt, (C.lb) 

where denotes a channel vector / for the fusion cen¬ 
ter that is randomly and independently generated at the 
(th iteration according to the known distribution of /, and 
(x,p; is the locally tight convex upper bound 

(IB. lb on the negative Bhattacharyya distance around the point 
Note that the cost function dC.lab depends on all 
the realizations of the channel vectors for I — 1,... ,j. 
The solution x*-®^ for the *th iteration of the outer loop is 
obtained by solving the problem (1C. lb iteratively over j, until 
a convergence criterion is satisfied. 

B. Optimization over p 

With the optimized waveform x = SSUM calculates 
the iterates by solving iteratively the following problems 

1 

argmin- ^7((')(x,p;/®|p(*’'-i)) (C.2a) 

P ^ 1=1 

s.t. 1 ^p<Pk, (C.2b) 

p„ > 0, n G Af, (C.2c) 


where ( x,p; ) is the convex upper bound ( IB. 3b 

on the negative Bhattacharyya distance around the point 

p(i.;-i) rpjjg iterate p^®^ is obtained by solving the problem 

(1C. 2b iteratively over j until convergence of the cost function. 

References 

[1] S. Jeong, O. Simeone, A. Haimovich, and J. Kang, “Optimization of 
multistatic cloud radar with multiple-access wireless backhaul,” in Proc. 
IEEE Radar Conf., Arlington, VA, May 2015, pp. 1650-1655. 

[2] M. Bemfeld, Radar signals: An introduction to theory and application. 
Elsevier, 2012. 

[3] A. W. Rihaczek, Principles of high-resolution radar. Wiley, 1967. 

[4] S. M. Kay, “Optimal signal design for detection of Gaussian point targets 
in stationary Gaussian clutter/reverberation,” IEEE Jour. Select. Topics 
in Sig. Proc., vol. 1, no. 1, pp. 31^1, Jun. 2007. 

[5] T. Kailath, “The divergence and Bhattacharyya distance measures in 
signal selection,” IEEE Trans. Comm., vol. 15, no. 1, pp. 52-60, Feb. 
1967. 

[6] S. M. Kay, “Waveform design for multistatic radar detection,” IEEE 
Trans. Aerosp. Electron. Syst., vol. 45, no. 3, pp. 1153-1166, Jul. 2009. 

[7] M. Naghsh, M. Modarres-Hashemi, S. Shahbazpanahi, M. Soltanalian, 
and P. Stoica, “Unified optimization framework of multi-static radar 
code design using information-theoretic criteria,” IEEE Trans. Sig. Proc., 
vol. 61, no. 21, pp. 5401-5416, Nov. 2013. 

[8] V. S. Chernyak, Fundamentals of multisite radar systems: multistatic 
radars and multistatic radar systems. Gordon and Breach Science 
Publishers, 1998. 

[9] M. Barkat and P. K. Varshney, “Decentralized CFAR signal detection,” 
IEEE Trans, on Aerospace and Electronic Systems, vol. 25, no. 2, pp. 
141-149. 

[10] R. Viswanathan and P. K. Varshney, “Distributed detection with multiple 
sensors: Part I - Fundamentals,” Proc. of the IEEE, vol. 85, no. 1, pp. 
54-63. 

[11] R. S. Blum, S. A. Kassam, and H. V. Poor, “Distributed detection with 
multiple sensors: Part II - Advanced topics,” Proc. of the IEEE, vol. 85, 
no. 1, pp. 64-79. 

[12] G. Alirezaei, M. Reyer, and R. Mathar, “Optimum power allocation in 
sensor networks for passive radar applications,” IEEE Trans. Wireless 
Comm., vol. 13, no. 6, pp. 3222-3231, Jun. 2014. 

[13] G. Alirezaei, O. Taghizadeh, and R. Mathar, “Optimum power allocation 
with sensitivity analysis for passive radai* applications,” IEEE Sensors 
Journal, vol. 14, no. 11, pp. 3800-3809, Jun. 2014. 

[14] China Mobile, “C-RAN: the road towards green RAN,” White Paper, 
ver. 2.5, China mobile Research Institute, Oct. 2011. 

[15] F. Gini and M. Rangaswamy, Knowledge based radar detection, tracking 
and classification. John Wiley & Sons, 2008, vol. 52. 

[16] J. R. Guerci, “Cognitive radar: A knowledge-aided fully adaptive ap¬ 
proach,” in Proc. IEEE Radar Conf, Washington, DC, 2010, pp. 1365- 
1370. 

[17] A. Gersho and R. Grey, Vector Quantization and Signal Compression. 
Boston: MA Kluwer, 1992. 

[18] T. M. Cover and J. A. Thomas, Element of Information Theory. John 
Wiley & Sons, 2006. 

[19] A. E. Gamal and Y.-H. Kim, Network Information Theory. Cambridge 
University Press, 2011. 

[20] R. Zamir and M. Feder, “On lattice quantization noise,” IEEE Trans. 
Info. Th., vol. 42, no. 4, pp. 1152-1159, Jun. 1996. 











13 


[21] M. W. Marcellin and T. R. Fischer, “Trellis coded quantization of 
memoryless and Gauss-Markov sources,” IEEE Trans. Comm., vol. 38, 
no. 1, pp. 82-93, Jan. 1996. 

[22] S. M. Kay, Fundamentals of Signal Processing-Estimation Theory. 
Englandwood Cliffs, NJ: Prentice Hall, 1993. 

[23] P.^.'W 2 i!s:?,hnQ:y, Distributed Detection and Data Fusion. Springer, 1997. 

[24] L F. Akyildiz, B. F. Lo, and R. Balakrishnan, “Cooperative spectrum 
sensing in cognitive radio networks: A survey,” Physical Communica¬ 
tion, vol. 4, no. 1, pp. 40-62. 

[25] S. P. Boyd and L. Vandenberghe, Convex optimization. Cambridge 
university press, 2004. 

[26] S. Khalili, O. Simeone, and A. M. Haimovich, “Cloud Radio-Multistatic 
Radar: Joint optimization of code vector and backhaul quantization,” 
IEEE Sig. Proc. Lett., vol. 22, no. 4, pp. 494^98, Oct. 2014. 

[27] D. R. Hunter and K. Lange, “A tutorial on MM algorithms,” The 
American Statistician, vol. 58, no. 1, pp. 30-37, Feb. 2004. 

[28] M. Razaviyayn, M. Hong, and Z. Luo, “A unified convergence analysis 
of block successive minimization methods for nonsmooth optimization,” 
SIAM Journal on Optimization, vol. 23, no. 2, pp. 1126-1153, Jun. 2013. 

[29] Z. Q. Luo, W. K. Ma, A. M. C. So, Y. Ye, and S. Zhang, “Semidefinite 
relaxation of quadratic optimization problems,” IEEE Signal Proc. 
Magazine, vol. 27, no. 3, pp. 20-34, May 2010. 

[30] M. Razaviyayn, M. Sanjabi, and Z.-Q. Luo, “A stochastic successive 
minimization method for nonsmooth nonconvex optimization with ap¬ 
plications to transceiver design in wireless communication networks,” 
arXiv preprint arXiv:1307.4457, Jul. 2013. 



